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Abstract. The spin-boson model has nontrivial quantum phase transitions at zero temperature induced by the spin- 
boson coupling. The bosonic numerical renormalization group (BNRG) study of the critical exponents /3 and S of this 
model is hampered by the effects of boson Hilbert space truncation. Here we analyze the mean-field spin boson model 
to figure out the scaling behavior of magnetization under the cutoff of boson states Nf,. We find that the truncation 
is a strong relevant operator with respect to the Gaussian fixed point in < s < 1/2 and incurs the deviation of the 
exponents from the classical values. The magnetization at zero bias near the critical point is described by a generalized 
homogeneous function (GHF) of two variables t = a - a^ and x = IjNi,. The universal function has a double-power 
form and the powers are obtained analytically as well as numerically. Similarly, m{a = a^) is found to be a GHF of e 
and X. In the regime .s- > 1/2, the truncation produces no effect. Implications of these findings to the BNRG study are 
discussed. 
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1 Introduction 

The spin-boson model (SBM) describing a quantum two-level 
system coupled to a dissipative environment appears in many 
areas of condensed matter physics ||T]|2l- It is a simple model 
for studying the dissipation and decoherence in a quantum sys- 
tem subjected to interactions with environment, such as a qubit 



The path integral with noninteracting blip approximation and 
adiabatic renormalization found no transitions in this case |2|. 
In contrast, the quantum- to-classical mapping theory for SBM 
showed that a quantum phase transition existed for Q < s < 
1 ||Tl|2]|8l and that the phase transition between the localized 
and delocalized phases should be in the mean-field universal- 
ity class in the regime < s < 1/2 [9,10|. Some important 



for the quantum computations |3|. The rich environment-induced properties of the equilibrium and non-equilibrium spin-boson 
quantum phase transitions in the spin-boson model also attract model has been obtained from the flow equation method ITTl 
much research attention recently. Experimentally, an environment lT2l[T3l . from which a first order phase transition was argued to 



induced localization transition has been observed in the Joseph- 
son junction systems Q, although typical solid-state two-level 
systems have a coupling strength below the critical threshold 
Q. Ref.|6| showed that optical forces allow realizing a vari- 
ety of spin-boson models, depending on the crystal geometry 
and the laser configuration. In most realistic situations the bo- 
son bath spectrum is super-Ohmic (s > 1) or Ohmic {s =1). 
Recently an experimental set up for realizing the spin-boson 
model in the sub-Ohmic {0 < s < I) and strong coupling 
regime is proposed in the mesoscopic metal ring systems Q. 
(s stands for the exponent of the bath spectral function). 

Many theoretical studies focus on the quantum phase tran- 
sitions in the spin-boson model [21 • For the super-Ohmic dissi- 
pation, the system is always delocalized and there is no phase 
transition. For the Ohmic dissipation, a quantum Kosterlitz- 
Thouless (KT) transition separates the delocalized phase at small 
dissipation from the localized phase at large dissipation, with 
critical coupling strength a » 1 in the small tunneling limit. 
The study of the SBM with sub-Ohmic bath is more difficult. 
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exist in the sub-Ohmic case. Recently, a perturbation approach 
based on unitary transformation was used to study the dynam- 
ics of SBM with sub-Ohmic bath and the results also suggested 
a first order localized-delocalized phase transition! 14 1. 

A comprehensive study on the critical behavior of SBM for 
the sub-Ohmic dissipation begins when the non-perturbative 
numerical renormalization group (NRG) method is extended 
to solve boson problems 1 1 5 , 1 6 , 1 7 1 . A continuous phase tran- 
sition is disclosed and the critical exponents obtained by NRG 
are non-classical and fulfill hyper-scaling relations. In the < 
s < 1/2 regime, the NRG results for the critical exponents 
/3 and 6 are different from the classical ones (J3 = 1/2 and 
6 - 3). This is in contrast to the claim of the quantum-to- 
classical mapping above ll9l [T0l . hence questioning the legiti- 
macy of the mapping theory in this model. Recently, the con- 
tinuous time Monte Carlo calculation was applied to the sub- 
Ohmic SBM |18|. It's found that in the regime < s < 1/2 
the critical exponents do take the classical values. A sparse 
polynomial space approach (SPSA) together with a standard 
exact diagonalization algorithm has been employed to calcu- 
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late the critical behaviors of the sub-Ohmic SBM lfT9l . The re- 
sults confirmed the Gaussian critical fixed point in the regime 
< s < 1/2. It is now realized that the exponents yS and 6 pro- 
duced by NRG are incorrect in the regime < s < 1/2, due to 
the boson state truncation used in the algorithm. How to avoid 
this error in the bosonic NRG is still an open question. 

The issue of quantum-to-classical mapping also appears in 
the Bose-Fermi-Kondo model (BFKM). In the spin-isotropic 
case, S. Kirchner et al. adopted a dynamical large-A^ limit of 
the SU(A^) method to study the effect of the Berry phase term of 
the spin path integral on the quantum critical properties. They 
attributed the emergence of the interacting fixed point to the 
interference of the Berry phase |20|. In their recent work, S. 
Kirchner et al. argued again that the presence of Berry phase 
changed the critical properties and claimed that the mapping 
theory failed for the sub-Ohmic BFKM |21 1. A scaling analy- 
sis for the Ising-BFKM was performed in Ref.[22], which was 
believed to have the same critical properties as the spin-boson 
model. The authors restated the failure of the quantum to clas- 
sical mapping for the QCP of the Ising-BFKM. They argued 
that the continuum limit taken in the mapping failed to capture 
the topological effect encoded in the Kondo spin-flips which 
was essential to the nature of the QCR 

Local boson Hilbert space truncation is frequently used in 
the numerical studies of boson systems, such as exact diagonal- 
ization (ED), NRG and ED-i-DMFT study of the Bose-Hubbard 
model ll23l . In these algorithms, calculating with a finite boson 
state number A^^ and then extrapolating the results to Nh - °° is 
believed to be sufficient to yield correct results. However, the 
NRG study of the spin-boson model in < 5 < 1 /2 regime 
shows that the simple extrapolation cannot guarantee the cor- 
rectness of the critical exponents p and cJ 1 24 1 . In order to figure 
out the role of A^^ in the critical behavior of the order param- 
eter m, in this paper, we present a numerical analysis to the 
mean-field spin-boson model. This model has a Gaussian criti- 
cal fixed point, the same as the spin-boson model in the regime 
Q < s < \ 12, and its exact critical exponents are known, i.e., 
/?= 1/2 and 5 = 3. 

The structure of this paper is as follows. In Section 2, we in- 
troduce the mean-field spin-boson model and its star-form. The 
methods of solution in the full as well as the truncated Hilbert 
space are presented. In Section 3, we present numerical results 
for the critical exponents /?, 6, and y. Scaling analysis of the 
magnetization function is employed to interpret these numer- 
ical results. The relations between our findings here and the 
NRG study of the spin-boson model are discussed. In Section 
4 we end with a brief summary. 



2 Model and Method 

2.1 The mean-field spin-boson hamiltonian 

The Hamiltonian of the spin-boson model reads {ti - I) |[l]|2l 



H = -jO-^ + -a-- + ^ Wjfljfl; + y ^ Ma\ + «,)• (1) 



between the two states. The environment is modeled as a col- 
lection of harmonic oscillators, which serve as the origin of 
dissipation! 1 2|. a', and a, are creation and annihilation opera- 
tors for the \-th phonon mode with frequency w,. /I, represents 
the coupling between the two-state system and the \-th phonon 
mode. The effect of the harmonic environment is characterized 
by the bath spectral function 



J{lj) - ^ / ^i S(Cl) - LJi), 



(2) 



which completely determines the influence of environment on 
the two level system. For simplicity, we use a power form of 
the spectral function with an energy cutoff u)c 



J{u) 



2naa) '^ cj^ 0{(jl>c 



■u). 



(3) 



Here Wc = 1 is the energy unit, or is a dimensionless coupling 
constant that characterizes the dissipation strength. The index 
s accounts for certain physical environment to which the two- 
state system couples. 

At zero temperature and zero bias, the competition between 
the quantum mechanical tunneling of the two states (leading to 
a delocalized phase) and the effect of the spin-bath coupling ( 
tending to localize the system into spin up or spin down state) 
leads to a phase transition between the two phases at a critical 
coupling ac (for a fixed A). 

We focus on the mean field spin-boson model. It is exactly 
solvable in the infinite boson Hilbert space, with an analytical 
expression of critical coupling Uc and the classical exponents. 
In the truncated boson Hilbert space, it can be solved numer- 
ically at high precision. It hence enables us to focus on the 
eff'ect of the local boson Hilbert space truncation on the criti- 
cal behavior. Carrying out the mean-field approximation to the 
spin-boson model, we obtain the mean-field Hamiltonian 






2 2 

-^2^,<« '"+«,.). (4) 



Neglecting the constant term, it can be written as the sum of 
decoupled Hamiltonians for isolated spin and displaced free 
bosons: 



^MF — ^snin "*" ^ boson-! 



with 



and 



'^snin — ni ^-V ' 



6 v~l Al -r- 

5 + Zj T<«'' + «'■> 



a-z. 



(5) 



(6) 



Here, the Pauli matrices a-f and cr, describe a two-state sys- 
tem, e is the energy difference and A is the tunneling strength 



Hbo 



2_^ w,a ! fl; + — ^ 2_j '^i ('^,' + «'■) ■ 



(7) 
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In the truncated boson Hilbert space, a, and a . are no longer 
canonical boson operators. One needs to resort to numerical 
calculations to solve Eq.(|5]l-(|7]l. Then, one has to specify the 
form of the SBM Hamiltonian, i.e., to parameterize the envi- 
ronment spectrum J(cij) and assign values to parameters /I, and 
(L>i. One common way to discretize and parameterize the bath 
degrees of freedom used in the NRG calculation is the loga- 
rithmic discretization. In order to make connection to the NRG 
studies, we use the same parametrization as in NRG. Follow- 
ing the procedure in Ref. (TT\ with an additional mean-field ap- 
proximation we arrive at the star-form mean-field Hamiltonian 
below 



with 



TJS, 



TTStar _ „ , 



_ TTStar . Tjstar 
~ spin boson ' 



9 +:r-py!r«<a« + «I) 

2 lyjJT^ 



(8) 



(r„ (9) 



and 



Here, the logarithmic discretization gives 



2na 



2 = ^[i_^-(i-)1a-''(i-V, (11) 

1 -I- 5 L J 



and 



1 + . 1-A-P"\ _„ 



(12) 



2 + s l-A-<i+-') 
A > 1 is the logarithmic discretization parameter. 

2.2 Numerical methods for truncated H'^p 

In the truncated boson Hilbert space, due to the decoupling of 
boson modes in Hf^^, it is possible to solve the boson part 
j^star Yiy exact diagonalization for each truncated mode. The 

boson •' ^ 

obtained boson average {a„ + a„) for n - 0, l,...,Nc are in- 
put into the spin Hamiltonian to solve for {cr,}. This iteration 
continues until convergence is reached. 

For the calculation of susceptibility x^ we start from the 
self-consistency equation of the order parameter m - {cr-,) at 
zero temperature. 



m(e,t,A) - 



zf2 + (e + r)2 +{e + t) V(e + t)^ + A^ 



-1, (13) 



with 



/ = t{e) 



V 



V^r 



J nan. 



(14) 



where 



l/h ^~ ^^ *^n I 



fli >, 



(15) 



and Nc boson modes are used. After some algebra, we obtain 
the final equation of the susceptibility 



X = -- 



(tl+A^ + t,^l^) ^ 



tt,+A^ 



A^ 



y-to ^tl 



i2 + 2f2 



^ 



1 W, N,, I ti, 2 

1 v-i v-i 2< gn\a„ + a„\k„ >^ 

^-nLLyn- 



=0 k=2 






(16) 



Here, fo is the tie) in Eq.(fT4l) at zero external bias. \g„ >, ef 
and \k„ >, e^ denote the eigenvectors and eigenvalues of the 
ground state and the A:,/, excited state of the n,;, boson mode, 
respectively. Due to the Hilbert space truncation, the summa- 
tion of k is limited to A^^ for each mode n. Eq.(fT6]l is evaluated 
numerically after the iterative solution is converged. 

We use total 101 different boson modes in the calculation, 
i.e., the summation of boson modes is cut off at Nc - 100. Due 
to the exponential decay of y„ and ^„, this summation is already 
numerically exact. For each boson mode, we retain Nh states. 
For simplicity, we use the boson number eigen states |0 >, |1 >, 
..., \Nh - 1 > as bases. In our calculation, the NRG parameter 
A = 2.0 is used. 



3 Results and Discussions 

3.1 Exact solution at A^i = oo 

In the full Hilbert space, a-, and a . are canonical boson opera- 
tors obeying the common commutation relation a,, a • - 5ij. 
One gets the self-consistent mean-field equations for Hmf 



and 



niF — e 
LA 



X^ - ^(mr - e)2 + /)2, 



with 






'0 T^(i> 



dix) — laojcl s. 



(17) 



(18) 



(19) 



This set of self-consistent equations can be solved and the crit- 
ical coupling strength for a given temperature T is 



sA 



n=0 



2tanh(/3zl/2)w, 
It reduces to ac - sAjiluc) at T = 0. 



(20) 
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To Study the quantum phase transition at T = 0, we focus 
on the following critical exponents y6, 6, and y, that are related 
to the behavior of the order parameter m lfT6ll . 



m{a > Uc, T - Q,e - Q) cc {a — Uc 



m(a - Qc, r = 0, e) oc |e 



IM 



Xiocia < ff,., r = 0, e = 0) oc (a^ - a) 



(21) 



(22) 



(23) 



The classical values /3 - 1/2 and 5 = 3 are obtained, as it 
should be. Via the partial derivative of m with respect to the 
local external field e, we get the zero temperature susceptibility 



X 



2ciJc{a-ac) ' 



4ciJc{ac-a) ' 



(a < Oc); 



(a > Q-c), 



(24) 



which gives y - I. 

Note that for the full boson Hilbert space, due to the linear 
nature of the logarithmic discretization and the transformation, 
the two forms of mean-field Hamiltonian, Hmf and HJ^^^ are 
essentially equivalent and belong to the same universality class. 
The mean-field approximation and the transformation can be 
interchanged in sequence. Therefore, the self-consistent mean- 
field equations Eq.([T7l[T8]l still hold for H 
placed by 



]^^'', but with r re- 



r 



rl 



V 2iL 



2QrWc(s + 2)[l-A-("+"]2 



(l + i)2(l-A-0[l-A-('+i>]' 
The zero temperature critical coupling Oc now reads 

"' ~ 2wc(s + 2)[l-(yV-('+i))]2 



(25) 



(26) 



In the limit A — > 1, o-c tends to the critical value adT = 0) = 
sA/{2(L)c), being consistent with Eq.(l20l). 



3.2 Analytical solution at Nh -2 
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Fig. 1. Exponent /? as a function of s at the parameters: A = 0.1, 
e = 0, A'i = 2. Dash line represents the analytical data; Olive squares 
symbol numerical data; Dash-dot line is for Ni, = oo. The inset shows 
the magnetization at i- = as a function of a with N^ = 2000. 



In the case of i = 0, we get a^ 
magnetization has the form: 



0. /3 diverges and the 



^ zJ/«A(l-i-A" ), 
^^ 8q'Wc(1-A-1) 



(28) 



The numerical results agree well with these results, as re- 
vealed in Fig. [1] The special case for 5 = is manifested as 
the inset in Fig.[T] In Fig. (l.a), the slope is approximately - 
0.02599 for A = 2.0, being consistent with the power of the 
exponential part in Eq.(l28ll. In Fig. (l.b), the slope is -0.5. 

Similarly, both analytical and numerical results for the crit- 
ical exponent 6 are available at A^^ = 2, as shown in Fig.|2l 



(s > k)\ 



rf' (0<.<i). 



(29) 



It is noted that in < i < 1/2, both/? and 6 for Nb -2 agree 
with the corresponding exponents in the spin-boson model ob- 
tained from NRG. A natural question is how the nonclassical 
critical exponents /3 and 6 in the regime < i < 1 /2 as given 
in Ea.(l27ll and Eq.(|29]l at Nh - 2 change to classical ones at 
Nb = oo. In the following section, we explore this issue at in- 
termediate Nb's. 



We obtained the critical exponents yS and 6 for H-^^^ analytically 
in the special case at A^^ - 2. (See in the Appendix for details.) 
We found that a^Nb = 2) = adNb = °°) as given in Eq.(l26]l. 
suggesting that ac is independent of the truncation Nb- This has 
indeed been observed in our numerical calculations at various 
Nh's- We also found for Nb -2 



y6 = 



2s ' 



(0 < . < i). 



(27) 



3.3 Numerical analysis at intermediate Nh's 

At intermediate A^^'s, we calculate < a„ + a„ > by exact di- 
agonalization and solve the mean-field equation iteratively. For 
s > 1/2, the numerical results always yield the classical crit- 
ical exponents /3 - 1/2 and 6 = 3, irrespective of the boson 
state truncation A^;/s. But for Q < s < 1/2, this is no longer the 
case. In Fig[3] we plot the dependence of the average magneti- 
zation m on the dissipation strength a for s = 0.2, at different 
A^i's {2 < Nb < 00). fff is calculated from Eq.(l26b at given 
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Table 1. Translation of m, t and x data lo Nt = 10 curve for s = 0.3. 




Fig. 2. Exponent 5"' as a function of s at the parameters: A = O.l, 
a = ttc, Nb = 2. The dash line is the analytical result; The olive squares 
are numerical data; The dash-dot line is for Ni, = oo. 
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15 


0.176 


0.20 


0.10 


0.57 


1.14 


0.50 


20 


0.301 


0.30 


0.15 


0.50 


1.00 


0.50 


30 


0.477 


0.48 


0.24 


0.50 


1.00 


0.50 


40 


0.602 


0.60 


0.29 


0.48 


1.00 


0.48 


80 


0.903 


0.86 


0.43 


0.48 


0.95 


0.50 



Here A is a positive number. Taking a logarithmic form of 
the equation above, we arrive at the following: 



m(T + Ot-A, X + a^A) — m(j, x) + a„,A, 



(31) 



where /=lg/, / - m, t, x, A. This equation implies that if m, t 
and X are shifted by Am, At, Ax, respectively, the curves will 
collapse. The ratios of the shifts give the critical exponents 
Am/AT — a,„/aj and Am/Ax — a,„/ax. Assigning A — —x/ox 
in Eq.dSTTi. we get 




Fig. 3. Igm v.s. Igr for different truncations at the parameters: A =0.1, 
s=0.2, e=0.0. The truncations are Ni,=2, 5, 10, 15, 20 and co from 
bottom to top. 



parameters. It is clearly seen that in the small a - a^ limit, m 
exhibits a perfect power law and the slope yS^ is identical to that 
of Nt - 2, namely ySrf = (1 - s)/{2s), deviating from the classi- 
cal exponent dramatically. Only for A^/, — > oo do we recover the 
anticipated mean-field result. In the upper region of the curve, 
the magnetization at finite but large A^^ tends to overlap with 
that at Nh - oo, giving a different slope jS„ = 1/2. As a result, 
two different power laws appear in the lower and upper regimes 
of the curve. In the following, we carry out a scaling analysis 
for m data with respect to parameters t = a - Oc and x = 1 /Nh 
on the basis of the generaUzed homogeneous function (setting 
6 = 0)123. 

We assume that the singular part of m is a GHF of t and x. 
That is. 



— / — ^T— „i — , ^ Am — 

m(T X, U) = m(T, X) X. 



(32) 



For simplicity, we denote z - t - —x. Fig. [3] suggests that 
we can assume a universal function as: 



m(T X, 0) — m(z, 0) — 



iPuZ + C„, (z » zo); 
[PdZ + Cd, (z-^zo). 



(33) 



Zo is some crossover value separating two regimes with differ- 
ent power laws. This gives 



m{T, x) 



Puf + I(f^ - f^A<) + C„, (r- » CX-); 
{PdT + I(f^ - ^Pd) + Cd, (t^ « cx^). 



(34) 



m{A"^T,A"'x) = A"-"m{T,x). 



(30) 



Cu, Cd and c are constants. The truncation independence of 

the upper power in Fig. [3] requires that the coefficient of x in 

_i_ _i_ 

T-T » cx"* regime should be zero, namely a,„/ax - f^uOr/ax 

or Pu — amlo-T- In the small t regime. Fig. [3] suggests that at 

finite truncations, fid is identical to that of Nh - 2, namely 

Pd^d- s)l{2s). 

Our assumptions in Eq.(l34l) are verified in Fig.|4|a). All the 
curves overlap with the curve at A^^ - 10 after proper trans- 
lations of m and 7. Translation details are in Tab. [T] We get 
Ami At - a„,/ar - 0.5, supporting yS„ = Pmf - 1/2, consis- 
tent with the specific case Nb = oo. According to Eq.(l34ll. the 
crossover point Tcr between the two power law regime yS„ = 1/2 
(upper power) andPd = (1 - s)/(2s) (lower power) is given by 
Tcr ~ x"^^"'. Tcr decHncs as Nh increases, leading to the expan- 
sion of the Pu = 1/2 regime. In the limit A^^, = oo, i.e., x = 0, 
Tcr is moved to zero, and the full mean field curve should be 
recovered. This is indeed observed in Fig. [3] 

In the regime r < Tcr, the second equation in Eq.(l34ll holds. 
Fig. Htb) shows Igm - yS^lgr v.s. Igx at fixed t in the lower 
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Fig. 4. (a) Overlap of curves with Nt = 15, 20, 30, 40, 80 to the curve 
with A'i = 10 after proper translation at .s = 0.3. (b) Scaling of the 
magnetization with respective to the truncation. From top to bottom 
s = 0.1,0.15,0.2,0.25,0.3, respectively. Here, we use lgT= -5.0 for 
i- = 0.1 and lgT=-6.0 for the others. Parameters arezl=0.1, 6=0.0. 



regime. Here /3d - (I - s)/(2s) is used. For various s values, 
the curves are linear with the same slope C. This confirms the 
second equation of Eq.(l34li and gives fidaj/a^^ - a„Jax - C, 
being independent of s. Numerical fitting gives the slope C ~ 
0.61 + 0.02. Taking into account of amldr - 1/2 and fid - 
(1 - s)/(2s), we get 






ICs 



_Cs_ 

l-2s' 



(35) 



(36) 



They are plotted as lines in Fig. |5] On the other hand, the 
average values of a,„/a^- and flr/flv obtained from translation 
of m data are plotted in Fig. |5]as functions of s. They match 
Eq.(l35ll and Eq.(l36ll remarkably well. 

Fig©-© show that even for a mean-field Hamiltonian, 
local boson state truncation can dramatically alter the criti- 
cal exponents, from the Gaussian exponents to the interacting 
ones. The boson state truncation, i.e., the artificial constraint 
of the local Hilbert space can be regarded as an additional lo- 
cal interaction introduced between bosons. For the mean-field 
model with Gaussian critical fixed point, this interaction be- 
comes dominant in low energies and leads the system to a new 
interacting critical fixed point. The critical exponent does not 
change continuously with the Hilbert space constraint. Instead, 
the constraint tunes the crossover point in the m-T curve below 
which the system flows into a truncation-dominated interacting 
fixed point. 

Our results shed some lights on the problem of bosonic 
NRG for the spin-boson model. ThatyS^ = (1 - s)/i2s) of H'^"; 
agrees with the /? obtained from NRG for the spin-boson model 
suggests that /3d itself is totally dominated by truncation. It is 
probable that the same scenario of A^^ also occurs in NRG. 



Fig. 5. Ratio of the scaling powers as functions of 5 at zl = 0. 1 . Orange 
triangles are a^/a^. = Am/Ax obtained from translation in Igm v.s. 
IgT figures at different s's as Fig. (|4]a) for s = 0.3; Dash line is the 
function a,„/a^. = -0.615/(2.? - 1). Wine squares are a^/a^ = AfjAx 
from translation in Igra v.s. Igr figures at different .s's like Fig. l4]a); 
The dot line is -1.22.s/(2.s- 1). Olive circles are aj/a, from translation 
in Igm v.s. Ige figures at different i's as Fig. ([6]b) for s = 0.2; The 
dash-dot line is -1.83.s/(2.s - 1). 



If this is true, it becomes clear why the bosonic NRG with 
boson state truncation cannot give correct classical exponent 
in this regime: the Gaussian critical fixed point in the regime 
< s < 1 /2 is overtaken by the strong relevant operator intro- 
duced by the Hilbert space truncation 11611241 . Moreover, one 
cannot improve the NRG critical exponent simply by increas- 
ing Nb if one focuses only on the small r limit. Instead, the cor- 
rect exponents can be crudely observed at finite energy scales 
from the truncated NRG calculation as shown in Ref.|24|. As 
demonstrated here, to fully disclose the role of Nb and to extract 
the accurate exponent/? from the bosonic NRG calculation, it is 
necessary to carry out a scaling analysis of Nb- It is worth men- 
tioning that the NRG study of the Ising-BFKM gives the inter- 
acting critical fixed point in the regime < s < 1/2, in contrast 
to the belief that flie Ising-BFKM and SBM belong to the same 
universality class II26II27I . In light of our study, we suggest that 
the NRG study for Ising-BFKM should be checked, using sim- 
ilar scaling analysis of Nb- 

In the following, we study the influence of truncation on 6- 
In Fig.|6] we show m as a function of bias e at the critical a^'s 
with different A^^'s and the collapse of these curves to the curve 
atNb = 10. Fig. |6la) shows that in the lowerregime the magne- 
tization has the same power dependence on e as that of Nb - 2, 
i.e., Sd = {I + s)/{l - s) as given in Eq.(|29]l. Fig.|6jb) shows 
the collapse of the Igm-lge curves under proper translations of 
X, e and m. This supports that, near the criticality m(£, x) is also 
aGHF: 



m{A"'e,A'''x) = i"'"m(e, jc). 
The ratios between the shifts are AmjAe 



(37) 



fl,„/fl(r and Am /Ax ■ 
Qin/cix- Details of the translation are in Tab.|2] According to Fig. 
|6ja) and Fig.l^b), m as a GHF of ;c and e can be formulated as: 
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Table 2. Translation of m, e and x data \.o Nt = 10 curve for s = 0.2. 



Nt 



Ax 



A1 



a„,lay a^ja^ o-mlo-t 



15 
20 
40 
80 
100 



0.176 
0.301 
0.602 
0.903 
1.000 



0.11 
0.19 
0.36 
0.56 
0.63 



0.036 

0.063 

0.12 

0.19 

0.21 



0.20 
0.21 
0.20 
0.21 
0.21 



0.63 
0.63 
0.60 
0.62 
0.63 



0.33 
0.33 
0.33 
0.34 
0.33 
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12 


-10 
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Fig. 6. Magnetization m as functions of bias e for different truncations 
at the parameters: A = id.\,a = a^, s = 0.2. (a) A',, = 2, 10, 40, 100, oo 
from bottom to top. (b) Overlap of curves with A'* = 1 5 , 20, 40, 80, 1 00 
to the curve with A'fc = 10 after proportional translation, (c) Scaling 
of the magnetization in terms of the truncation, i- = 0. 1, 0. 15, 0.2, 0.3 
from top to bottom, y = [m - 6(1 - i')/(l -I- s)\l\lCsl{\ + 5)], here we 
use ? = -9.0 and C = 0.61. 



m(e, x) — 



6^,^e + Cu, 



(e°e » cx-j ); 



(38) 



lV^ + ^(f:-|V) + C,, {e-«cx-). 



We get 5^' = flm/fle - Ami At ^ 1/3, signaling a clas- 
sical power law in the large e regime. In Fig. |6jc) we plot 
ijn - dj^Tj/iam/ax - 5^'fle/flx) versus x for different s. They 
are linear functions of x with an average slope 1.01 ± 0.05, as 
expected from the second equation of Eq.(l38ll. 

In Fig|5] flf/flj: obtained from the translation procedure is 
plotted as circles, to be compared with the analytical curve 
cie/cix = 3Ci/(l -2s). Here we use the value C - 0.61 obtained 
previously. They agree very well. These numerical results con- 
firm our assumption in Eq.(l38ll. 

This shows that m(e, x) is also a function with two different 
power law regimes: the low energy regime with l/Sd and the 
upper one with 1 /6mf = 1/3. They are separated by a crossover 
point ecr ~ x"^!"', similar to the situation in mij, x). 

In the following, we discuss the susceptibility exponent y. 
In Fig. I2I 1 /a' versus a is plotted. The zero point in the figure is 
the critical point. Within the critical region, the susceptibility 
exponent remains the classical one, y - I under finite boson 
state truncation in < 5 < 1/2. This shows that the boson state 



0.000 




-0.002 



0.0102 



0.0103 



0.0104 



0.0105 



a 



Fig. 7. Inverse of susceptibility l/x ^^ a function of a for different 
truncations Ni, at the parameters: A = 0.1, s=0.2, e = 0. The zero 
point of 1/x gives the critical coupling strength a^. Nt = 2, 5, 10, 15, 
20, 40, 80, 100, 00 from top to bottom, respectively. 



truncation does not influence y. Similar observations are made 
in the NRG, MC, SPSA and extended coherent state approach 
studies LjJIUlIigill. 

From the analysis above, we observe that in the low energy 
regimes, i.e., t'^"' «c cx'^"' and e'^"' <K cx'^"', fi and 6 are no 
longer determined by the scaling exponents Um, Ot, Oe- Instead 
new parameters are needed to define those critical exponents. 
We also find that in both the high and low energy regimes, the 
scaling relation /3(i5 - 1) = y is fulfilled. 

The same analysis is caiTied out at a different A and A: 
A - 0.2 and A - 3.0. Conclusions are qualitatively the same. 
From the analysis above, we believe that it is not the failure 
of quantum-to-classical mapping for the spin-boson model, but 
the numerical concomitant, i.e., the truncation that leads to new 
exponents in the regime < s < 1/2. Our work shows that the 
boson state truncation Nt severely changes the critical behavior 
of the mean-field spin-boson model: it plays the role of a new 
scaling parameter and changes the low energy exponent/? and 
6 into a truncation-dominated exponent (Eg.dZTll and Eq.(|29ll). 



4 Summary 

To conclude, we have analyzed the mean-field spin boson model 
under boson state truncation. We focus on the quantum critical 
behavior of the model and caiTy out the scaling analysis of the 
magnetization in terms of truncation parameter Nh- Our work 
shows that for the mean field spin boson model the truncation 
gives rise to a relevant perturbation in the regime < s < 1/2. 
The effect of the truncation, while cutting no ice for the expo- 
nents /3 and 6 in the regime s > 1/2, dominates in < .? < 1/2. 
As is shown in in Fig.|3]and in Fig.|6] the magnetization m un- 
der truncation is a generalized homogeneous function in the 
regime < s < 1/2: in the high energy regime, m fulfills 
a classical power law; in the low energy regime it exhibits a 
double power behavior, and the truncation A^^, becomes a new 
scaling parameter, besides the coupling strength t and exter- 
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nal field e (Eq.(l34b and Eq.(l38ll). Moreover, we find that as A^^ 
increases, the crossover point between two power law regimes 
moves downwards and to zero energy scale at A^^ = oo (namely 
without boson states truncation). Implications of these findings 
to the bosonic NRG study are discussed. For the spin-boson 
model within < i < 1 /2, it is now believed that the low en- 
ergy critical exponents yS = (1 - s)/{2s) and S = {I + s)/{l - s) 
obtained by NRG are also due to boson state truncation A^^ [ 16; 
l24ll . We conjecture that similar scaling behavior of A^^, should 
also occur in the NRG study of the spin-boson model. A cor- 
rect and accurate extraction of/3 and 6 for the spin-boson model 
from NRG requires similar scaling analysis as done here. Fur- 
ther studies in this direction are in progress. 

N. H. Tong acknowledges helpful discussions with Ralf 
Bulla and Matthias Vojta. This work is supported by National 
Basic Research Program of China (grant number 2007CB925004), 
and by the NSFC under grant number 10674178. 



A Appendix 

In the Appendix we give the analytical derivation of the critical 
exponents at Nb - 2. Following similar procedures in Sec. 2.2, 
we can get the self-consistent equations for the magnetization 
m at A^i = 2 



^'^-^tZ'^-^ z 



d„]. 



«— HO + l 



When n <K no, i.e., ^^ » ^^, 






2 2 

my; 



CfS 



2k 



k=0 



On the contrary, when n » riQ, i.e., ^„ <K — — 



z)2 + (e + e' -i-f)2' 



(39) 



dn 



^Jn 



(1 + 



f«^ ^-I 



m^yl 



Z n[lk{ 

QkA- 



k=0 



= ^fsTT and Qk = 

.[_1 _(fc_l)](C-^2 



-2A:-1 



c:"V+"2b2' 



(45) 



(46) 



with C] 



(47) 



-1/2 



Here, Pk 

_i(_i _ 1) 

B = ^„yV". Replacing J„ in Eq.lgsl) with Eq.gSll and Eq.gTji, 

respectively, we arrive at the summation: 



Xit-i , "'"■ ^k 
= 1), a2 = y2A"('+''> and 



f = V(e + e')2-Hzf2, 



(40) 



m \-~t , 



m^'' -Tm-^]. 



(48) 



e = ) d„. 



dn = 



r,^ 



Vi^ 



(41) 



(42) 



with y„ and ^„ satisfying Eq. (fTTl i and Eq. (fT2l l. respectively. 
Since we are interested in the critical exponent/?, first we take 
the bias e to be zero. Combining Eq. (l39b and Eq.(l40ll. we get 
the following equation: 



Here, f^ = P^/[l - A*^*^' "'*"*]. It is hard to get the analytical 
form for T, anyhow, it does not contribute to the exponent of 
the order parameter. Neglecting the trivial solution m - 0, we 
get the following self-consistent equation by substituting e' in 
Eg.(l43ll with Eq.(l48b. 



[^ P[m^'' - TmT^'fil - m^) = n^A^ 



(49) 



Considering that the order parameter approaches zero near 
the critical coupling a^ , we get 



P[^ + P\m^ - Tm'^' =nA. 



(50) 



V? 



+ /12 



(43) 



The final critical coupling strength Uc is 



The newly obtained Eq.(l43Tl together with Eq.flTt com- 
poses the self-consistent equations. To calculate the summation 
in Eq. lHTt . we define a parameter no determined by the follow- 
ing equation: 






2 2 



(44) 



Then Eq. (l4TI) can be divided into two parts: 



zl(.;-H 1)2(1 -A-')[l-(A-^'+2))] 
2w,(iH-2)[l-(A-<'+i>)]2 ' 



(51) 



the same as in Eq.(l26]l for Nb - c». For s > 1 12, the third 
term in Eq.(l50t is negligible and one gets m oc (a - a'e)'/^- for 
s = 1/2, the second and third term have the same power, i.e., 
(1 - s)/{2s) = 1/2 (s=l/2), so one gets m cc (a - ffc)'^^; for 
< 5 < 1/2, ignoring the second term in Eq.dSOll. one gets 
m cc (a - a^)*' ■'^^*^*'- In summary, for A^;, = 2, we get 
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I 

2' 



1-.; 

2.! ' 



(s > i); 



(0 < 5 < i). 



(52) 



In the case of s = 0, a small difference lies in Eq.(l46]l at 
k = 0. If s = Q and k = Q, the common ratio of the geometric 
series in Eq.(l45b for n < no, namely /iW"^)"'l, is unity and 
the general summation formula is not applicable. Taking this 
speciality into account and following similar procedure of the 
case i ?i 0, one gets the critical behavior of the magnetization 
as 
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J AlnAj 1 +A ) 

in oc a~^e sn-wcd-zi-') _ 
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As far as 6 is concerned, similar analysis gives the follow- 
ing expression 



5 = 



1-.S' 



is > i); 



(0 < . < i). 



(54) 



References 



1. U. Weiss, Quantum Dissipative Systems (World Scientific, Singa- 
pore 1993). 

2. A. J. Leggett, S. Chakravarty, A. T Dorsey, M. P. A. Fisher, A. 
Garg and W. Zwerger, Rev. Mod. Phys. 59, 1 (1987). 

3. F. Verstraete, D. Porras and 1. 1. Cirac, Phys. Rev. Lett. 93, 227205 
(2004). 

4. I. S. Pentilla, U. Parks, P. I. Hakonen, M. A. Paalanen, and E. B. 
Sonin, Phys. Rev. Lett. 82, 1004 (1999). 

5. T Hayashi, T Fujisawa, H.D . Cheng, Y H. Jeong, and Y Hi- 
rayama, Phys. Rev. Lett. 91, 226804 (2003). 

6. D. Porras, F. Marquardt, I. von Delft, and I.I . Cirac, Phys. Rev. A. 
78, OIOIOI(R) (2008). 

7. N. H. Tong and M. Vojta, Phys. Rev. Lett. 97, 016802 (2006). 

8. F. I. Dyson, Commun. Math. Phys. 12, 91 (1969). 

9. M. E. Fisher et al., Phys. Rev. Lett. 29, 917 (1972) 

10. E. Luijten and H. W. I. Blote, Phys. Rev. B. 56, 8945 (1997) 

11. S. K. Kehrein and A. Mielke, Phys. Lett. A. 219, 313 (1996). 

12. T Stauber and A. Mielke, Phys. Lett. A. 305, 275 (2002). 

13. T Stauber, Phys. Rev. B. 68, 125102 (2003). 

14. Z. Lu and H. Zheng, Phys. Rev. B. 75, 054302 (2007). 

15. R. Bulla, N. H. Tong and M. Vojta, Phys. Rev. Lett. 91, 170601 
(2003). 

16. M. Vojta, N. H. Tong and R. Bulla, Phys. Rev. Lett. 94, 070604 
(2005). 

17. R. Bulla, H. I. Lee, N. H. Tong, and M. Vojta, Phys. Rev. B. 71, 
045122 (2005). 

18. A. Winter, H. Rieger, M. Vojta and R. Bulla, Phys. Rev. Lett. 102, 
030601 (2009). 

19. A. Alvermann and H. Fehske, Phys. Rev. Lett. 102, 150601 
(2009). 

20. S. Kirchner and Q. Si, larXiv:0808.2647V l icond-mat.str-ell 
(2008). 

21. S. Kirchner and Q. Si , Physica B, 404 2904 (2009). 



